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Abstract 

Freshwater fish evolution has been shaped by changes in the earth's surface 
involving changes in the courses of rivers and fluctuations in sea level. The 
main objective of this study is to improve our knowledge of the evolution of 
loricariids, a numerous and adaptive group of freshwater catfish species, and 
the role of geological changes in their evolution. We use a number of different 
phylogenetic methods to test the relationships among 52 representative taxa 
within the Neoplecostominae using 4676 bps of mitochondrial and nuclear 
DNA. Our analysis revealed that the subfamily Neoplecostominae is monophy- 
letic, including Pseudotocinclus, with three lineages recognized. The first lineage 
is composed of part of Pareiorhina rudolphi, P. cf. rudolphi, and Pseudotocinclus; 
the second is composed of Isbrueckerichthys, Pareiorhaphis, Kronichthys, and the 
species Neoplecostomus ribeirensis; and the third is composed of Pareiorhina car- 
rancas, P. cf. carrancas, Pareiorhina sp. 1, a new genus, and all the species of 
the genus Neoplecostomus, except N. ribeirensis. The relaxed molecular clock cal- 
ibration provides a temporal framework for the evolution of the group, which 
we use for a likelihood-based historical biogeographic analysis to test relevant 
hypotheses on the formation of southeast Brazil. We hypothesize that headwa- 
ter capture events and marine regressions have shaped the patterns of distribu- 
tion within the subfamily Neoplecostominae throughout the distinct basins of 
southeast Brazil. 



Introduction 

The rivers of the Neotropics are hugely diverse, estimated 
to contain more than 7000 fish species (Eschmeyer 2011). 
Interestingly, most Neotropical ichthyofauna diversified 
within a comparatively short timeframe (post K/T bound- 
ary) (Albert and Reis 2011). The drainages of southeast 
Brazil, comprising the Upper Rio Parana Basin and coastal 
rivers (Rio Paraiba do Sul and Rio Ribeira de Iguape), are 
among the most species-rich regions with regard to 
freshwater fish (Abell et al. 2008). Among Neotropical 
freshwater fishes lineages, the Loricariidae represent the 
largest family with 973 species recognized (Eschmeyer 
2011) and at least 300 waiting to be described (Reis et al. 
2003). Within Loricariidae, Neoplecostominae include 37 
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comparatively small-bodied species (Eschmeyer 2011) in 
the genera Kronichthys, Isbrueckerichthys, Pareiorhaphis, 
Pareiorhina, Neoplecostomus, and Pseudotocinclus (Montoya- 
Burgos et al. 1998; Armbruster 2004; Chiachio et al. 
2008). Neoplecostominae representatives are restricted to 
the Southern and Southeastern regions of Brazil, making 
them an ideal group to test hypotheses of historical bio- 
geography relevant to these regions. These fish are found 
in small-to-medium-sized streams with clear shallow 
water (commonly <1 m deep), moderate-to-strong cur- 
rents, on loose stones, and sandy substrates (Langeani 
1990; Pereira and Reis 2002). The Neoplecostominae have 
a long and complex history of taxonomy and systematics, 
with morphological and molecular studies focusing on 
evolution of the subfamily and also in Loricariidae as a 
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whole (Eigenmann and Eigenmann 1890; Regan 1904; 
Gosline 1947; Isbrucker 1980; Howes 1983; Schaefer 1987; 
Montoya-Burgos et al. 1998; Armbruster 2004; Pereira 
2005; Reis et al. 2006; Chiachio et al. 2008; Cramer et al. 
2008, 2011). However, the historical evolutionary rela- 
tionships among endemic species from the southeast of 
Brazil and the mechanisms by which they diversified in 
space and time remain understudied. This has been a chal- 
lenge for generations of evolutionary biologists, mainly 
because of the difficulties in reconstructing the complex hy- 
drogeological processes that have shaped the Neotropics 
(Ribeiro et al. 2006; Albert and Reis 2011). 

Despite of the complexity of biotic spatiotemporal dis- 
tribution, in the last few decades, our understanding of 
the historical biogeography of Neotropical freshwater 
fishes has improved with the advent of plethora of 
phylogenetic, paleontological, paleoclimatic, and paleoen- 
vironmental data (Lundberg and Chernoff 1992; Casciotta 
and Arratia 1993; Gayet 2001; Gayet and Meunier 2003; 
Lundberg and Aguilera 2003; Lundberg 2005; Sanchez- 
Villagra and Aguilera 2006; Malabarba and Lundberg 
2007; Malabarba and Malabarba 2008; Albert and Reis 
2011). Abell et al. (2008), in an extensive synthesis on 
freshwater fish endemism throughout the world, proposed 
the following four eco-regions for southeast Brazil: the 
coastal drainages of Paraiba do Sul, Ribeira de Iguape and 
Southeastern Mata Atlantica, and the interior drainage of 
Sao Paulo State, known as Upper Parana. Such eco- 
regions are supported by published phylogenetic data 
indicating distinct ichthyofaunal divisions between the lit- 
toral drainages and Upper Parana (Montoya-Burgos 2003; 
Hubert and Renno 2006). In a recent study, Chiachio 
et al. (2008) tested biogeographic hypotheses involving 
the subfamilies Neoplecostominae, Otothyrinae, and 
Hypoptopomatinae using two dispersal-extinction-clado- 
genesis models to infer the spatiotemporal distribution of 
these fishes. Maximum-likelihood (ML) reconstructions 
of ancestral ranges indicated a marked division between 
the Amazonian origin of the Hypoptopomatinae and the 
eastern coastal Brazil + Upper Parana origin of 
the Neoplecostominae and Otothyrinae. 

Ribeiro (2006) has provided an overview of the princi- 
pal aspects of geological history along the continental 
margin of eastern South America and explained the 
biogeographic history of fish fauna from the Brazilian 
crystalline shield, focusing mainly on the coastal drainages 
of eastern Brazil. According to Ribeiro (2006), many clad- 
ogenetic events associated with tectonics and erosive 
processes (which are still active today across eastern South 
America) may be influencing ichthyofaunal distribution 
and speciation (Fig. 1). Ribeiro (2006) described river 
headwater capture events as a complex process associated 
with the main fault systems, which are more susceptible 



to erosion. In this process, portions of the tributaries of a 
river in one hydrographic basin are "captured" by adja- 
cent basins resulting in isolated populations (Fig. 1). This 
process occurred several times in the formation of South- 
eastern Brazilian shield. Albert and Reis (2011) suggested 
that between 15 and 28 million years ago several head 
water captures occurred in the Upper Parana basin and 
coastal rivers in Tremembe formation, whereas Ribeiro 
(2006) suggested the same for the Ponta Grossa Arch 
formation. Also, Montoya-Burgos (2003) studying the 
biogeography of the loricariid genus Hypostomus sug- 
gested that headwater captures occurred between the 
Upper Rio Parana and Eastern Coastal Rivers 4.2 million 
years ago. These vicariant processes were followed by allo- 
patric divergence and speciation within Hypostomus. A 
similar hypothesis has been proposed by Weitzman et al. 
(1988) working with Glandulocaudini species. 

Several authors have suggested that species diversity in 
Brazil could also be related to the occurrence of sea level 
fluctuations that have allowed the dispersal between 
neighboring rivers via delta connections or coastal 
marshes (Hoorn 1993; Monsch 1998; Hernandez et al. 
2005; Roddaz et al. 2005; Rebata et al. 2006; Hoorn and 
Wesselingh 2010). Montoya-Burgos (2003) used this class 
of events to explain the split between the lineages inhab- 
iting the eastern coastal rivers and those inhabiting 
Amazon drainages for species of the genus Hypostomus 
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Figure 1. (a) Scheme activation of ancient faults and the 
erosive process in Serra do Mar formation (arrow indicate the 
direction of the flow); (b) Scheme of the headwater captures resulted 
of the activation ancient faults and erosive process. Modified from 
Almeida and Carneiro (1998) and Albert and Reis (2011). 
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between 5 and 6 million years ago. Bloom and Lovejoy 
(2011) further suggested that vast areas of South Amer- 
ica, particularly, the Amazon basin, have relatively low 
elevations (<100 m) and are thus likely to be impacted 
by sea level fluctuations. Moreover, from the middle of 
the Miocene to the present, the climate has become 
colder and has undergone periodic climatic oscillations 
that were responsible for marine transgression and 
regression. Albert and Reis (2011) suggest that large mar- 
ine regressions occurred between 76-80, 66-71, 45-50, 
35-42, 23-34, 9-12, and from 3.5 million years ago to 
the present when climatic oscillations intensified. The 
main factors responsible for climatic oscillation are 
atmospheric C0 2 concentration, thermohaline circula- 
tion, the positions of the continental landmasses, and the 
tilt of the planet's axis (Budyko 1969; Petit et al. 1999; 
Rahmstorf 2003; Rothman et al. 2003). 

In this study, we use a phylogenetic approach to recon- 
struct a species tree for the Neoplecostominae and infer 
their age, origin, and ancestral dispersal routes through- 
out southeastern Brazil. We test relevant biogeographic 
hypotheses (Table 1) using a dispersal-extinction- 
cladogenesis approach, and discuss patterns likely to have 
influenced the spatiotemporal distribution of Neotropical 
fish fauna in the Southeastern part of Brazil. 

Material and Methods 
Sampling 

All fishes collected for this study (Table SI) were collected 
in accordance with Brazilian laws, under a permanent 
scientific collection license in the name of Dr. Claudio 
Oliveira. Additionally, our laboratory has special federal 



permission to keep animals and tissues from a public 
collection under our care. To work with the animals, we 
follow all the ethical prescriptions stated by our internal 
committee of ethic involving animal experiments (proto- 
col number 388). The species were collected using hand 
nets, from a variety of locations across southeastern Brazil 
(Fig. 2, Table SI). After collection, the animals, to do not 
suffer any kind of pain, were anesthetized using 1% 
benzocaine in water and either preserved in 95% ethanol 
for molecular studies or fixed in 10% formaldehyde for 
morphological studies. Vouchers of all sequenced samples 
were deposited in the collection of either the Laboratorio 
de Biologia e Genetica de Peixes (LBP), Departamento de 
Morfologia, Instituto de Biociencias, Universidade Estadual 
Paulista, Botucatu, Sao Paulo, Brazil, the Nucleo de 
Pesquisas em Limnologia, Ictiologia e Aquicultura (NUP), 
Universidade Estadual de Maringa, Parana, Brazil, or the 
Museum of Natural History of the City of Geneva 
(MHNG), Geneva, Switzerland. In the course of this 
study, several new species and two new genera were 
found and will be further described in subsequent 
publications. 

DNA sequencing 

Total DNA was extracted from ethanol-preserved muscle, 
fin, and liver samples using the Wizard Genomic DNA 
Purification Kit (Promega, Madison, Wisconsin). Partial 
sequences of 16S rRNA (700 bp), 12S rRNA (900 bp), 
cytochrome c oxidase subunit I (COI, 700 bp), cytochrome 
b (Cytb, 900 bp), and F-reticulon 4 (1900 bp) were ampli- 
fied using polymerase chain reaction (PCR) with the primers 
described in Table S2. The primer concentration was 
5 pmol/^L. Mitochondrial genes were amplified with a 



Table 1. Events which have influenced diversification of fish species in Southeastern Brazil. 



Event 



Predictions 



Methods used 



References 



Boundary displacements between 
the Upper Parana and the Littoral 
Drainages 

Head water capture events 
between major drainages 



3. Miocene, Pliocene and Pleistocene 
climate oscillation 



Sister lineages from different 
endemic drainages will be 
monophyletic 

(i) Clades (intraspecific lineages or 
species groups) will not be 
monophyletic for a given drainage; 
(ii) lineages in a given 
watershed will be nested within 
the phylogeny of a larger clade 
located in a different watershed 
with contiguous headwater. 

(i) Delta connection in coastal rivers 
in marines regression possibility 
the dispersal of species to others 
drainages. 



Phylogeny 



Phylogeny and molecular 
clock 



Phylogeny and Molecular 
clock 



Present study 



Weitzman et al. (1988), Ribeiro 
et al. (2006), Lundberg (1998), 
Montoya-Burgos (2003), Hubert 
and Renno (2006) 



Montoya-Burgos (2003), Hoorn and 
Wesselingh (2010), Abell et al. 
(2008), Hubert and Renno (2006), 
Ab'Saber 1979), Colinvaux et al. 
(2000), Petit et al. 1999), 
Wilkinson et al. (2006) 
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total volume of 25 fiL for 35 cycles (30 sec at 95°C, 
45 sec at 48-54°C, and 80 sec at 72°C). The nuclear 
markers were amplified in two PCR experiments; the first 
amplification using the primers Freticul4-D and Freticul4-R 
with a total volume of 12.5 for 37-40 cycles (30 sec at 
95°C, 30 sec at 48°C, and 135 sec at 72°C); and the 
second amplification using the primers Freticul4 D2, Fre- 
ticul4 R2, and Freticul4 iR with a total volume of 12.5 [iL 
for 37-40 cycles (30 sec at 95°C, 30 sec at 53-54°C, and 
135 sec at 72°C). The PCR was amplified using the Gotaq 
commercial kit (Promega). After that, the products were 
identified on a 1% agarose gel. The PCR products 
were purified using ExoSap-IT® (USB, Affymetrix Cor- 
poration, Cleveland, Ohio) following the manufacturer's 
instructions. The purified PCR products were used to 
make a sequencing PCR using the BigDye M Terminator 
v 3.1 Cycle Sequencing Ready Reaction Kit (Applied 
Biosystems- Life Technologies do Brasil Ltda, Vila Guar- 
ani, SP, Brazil). Subsequently, the amplified DNA was 
purified again and loaded on a 3130-Genetic Analyzer auto- 
matic sequencer (Applied Biosystems). All sequences obtained 
in this study were deposited in the GenBank (Table SI). 

Phylogenetic analyses 

An Astroblepus sp. (Astroblepidae) was used to root the 
Neoplecostominae phylogeny. Additionally, samples of 
Delturinae (Hemipsilichthys gobio and H. papillatus), 
Loricariinae (Rineloricaria jaraguensis), Hypoptopomatinae 
(Hypoptopoma inexpectatum), Otothyrinae (Corumbataia 
cuestae), and Hypostominae (Hypostomus nigromaculatus) 
were included in the analysis as additional outgroups. 
Sequences from databases was not included in our analy- 
sis because only sequences of the gene F-reticulon 4 are 
available, and the increase in gaps number should reduce 
the phylogenies accuracy like suggested by Dwivedi and 
Gadagkar (2009). Also, we have a great sample of 
Neoplecostominae species and that was considered mono- 
phyletic by Chiachio et al. (2008) using one of the molec- 
ular markers included in our analysis, F-reticulon 4. All 
individual sequences for each species were initially ana- 
lyzed using the software BioEdit 5.0.9 (Hall 1999) and a 
consensus sequences were obtained. Afterward, all 
sequences were independently aligned using the software 
Muscle (Edgar 2004). Nucleotide variation, substitution 
patterns, and genetic distances were examined using 
MEGA 5.0 (Tamura et al. 2007). All different nucleotide 
substitution models were tested using Modeltest 3.7 
(Posada and Crandall 1998) to 19 different partitions of 
the matrix. For each phylogenetic program analysis, a 
related model was used when the program do not present the 
nucleotide substitution models found for Modeltest. The parti- 
tion and the models utilized for each one are shown in Table S3. 
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Maximum-likelihood analyses were performed using 
RAxML Web-Servers (Stamatakis et al. 2008). It imple- 
ments a faster algorithmic of heuristic search with 
bootstrap pseudoreplicates (RBS). Bootstrap resampling 
(Felsenstein 1985) was applied to assess support for 
individual nodes using 1000 replicates. Random starting 
trees were used for each independent ML tree search and 
all other parameters were set on default values. All ML 
analyses were conducted under GTR, as RAxML 
only applies this model, using or not G according to 
Modeltest 3.7 (Posada and Crandall 1998) search 
(Table S3). 

Bayesian inference (BI) (Huelsenbeck and Ronquist 
2001) was performed evaluating alternative tree topolo- 
gies through the estimation of posterior probabilities 
using MrBayes v.3.0 (Ronquist and Huelsenbeck 2003). 
MrBayes implements three main types of nucleotide sub- 
stitution models and they are selected by setting the num- 
ber of substitution types using lset nst to 1, 2, or 6. This 
parameter was selected to each matrix partition 
(Table S3). Four chains were run simultaneously for 
50,000,000 generations and every 1000th generation, a 
tree was sampled. The above analysis was performed 
twice. The distribution of log likelihood scores was exam- 
ined to determine stationarity for each search and to 
decide if extra runs were required to achieve convergence, 
using the program Tracer 1.4 (Rambaut and Drummond 
2007a). All sampled topologies beneath the asymptote 
(5,000,000 generations) were discarded as part of a burn- 
in procedure, and the remaining trees were used to con- 
struct a 50% majority-rule consensus tree in Paup* 
(Swofford 2003). 

Maximum-parsimony (MP) was performed using the 
software TNT (Goloboff et al. 2008). No a priori 
weighting or ordering of character states was used and 
gaps were treated as missing data. Phylogenies were 
constructed under the "new technology search" meth- 
odology (Goloboff 1996, 1999), using the options 
"sectorial search", "ratched", "drift", and "tree fusing" 
with their default values and employing a driven search 
with initial level setting at level 100 and checking level 
at every two hits. Consistency and retention indexes 
(RIs) were calculated with the script "stats" of TNT. 
Clade robustness was assessed using 1000 bootstrap (B) 
pseudoreplicates (Felsenstein 1985) with the same 
parameters cited above. Bremer support values (BS) 
(Bremer 1988) were calculated with the script "bremer" 
script for TNT. 

Molecular clock and historical biogeography 

The uncorrelated relaxed molecular clock (lognormal) 
was calibrated using BEAST (Bayesian evolutionary analy- 
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sis sampling trees) vl.6.2, (Drummond and Rambaut 
2007) on the same dataset used above with the same par- 
titioning scheme (Table S3). The models utilized for each 
partition are shown in Table S3. We used the origin of 
Loricariidae family (85-90 million years ago), as esti- 
mated by Lundberg et al. (2007), for our calibration point 
at the root of the tree. This was implemented with a log- 
normal prior offset of 85 million years ago with a mean 
and standard deviation of 1. We used a Birth-Death 
model for speciation likelihood and a starting tree 
obtained from ML analysis. The analysis was run for 100 
million generations and sampled every 1000th generation. 
Stationarity and sufficient mixing of parameters 
(ESS > 200) was checked using Tracer vl.5 (Rambaut and 
Drummond 2007a). A consensus tree was built using 
TreeAnnotator vl.6.2 (Rambaut and Drummond 2007b). 

Maximum-likelihood inference of geographic range 
evolution was performed using the Dispersal-Extinction- 
Cladogenesis (DEC) as implemented in Lagrange v 2.0 
(Ree et al. 2005; Ree and Smith 2008). The DEC model 
specifies instantaneous transition rates between discrete 
distribution ranges along the branches of a phylogenetic 
tree and uses it to assess likelihoods of ancestral distri- 
bution ranges at cladogenesis events (Ree et al. 2005; 
Ree and Smith 2008). We used the unconstrained DEC 
model because the geographic region of interest was 
only divided into two areas (interior and littoral drain- 
ages). We classified the divisions of endemic areas 
according to the eco-region classifications of Abell et al. 
(2008), Montoya-Burgos (2003), and Hubert and Ren- 
no (2006). "Littoral drainages" was a classification 
composed of the following eco-regions: Northeastern 
Mata Atlantica, Paraiba do Sul, Ribeira do Iguape, 
Southeastern Mata Atlantica, and Fluminense. "Interior 
drainages" was composed of Upper Parana, Iguassu, 
and Sao Francisco. 

Results 

Phylogeny of the Neoplecostominae 

Partial sequences of four mitochondrial genes (12S rRNA, 
16SrRNA, COI, Cytb) and one nuclear gene (F-reticulon 
4) were obtained from 52 specimens representing 47 
loricariid species and one Astroblepidae species 
(Table SI). The combined sequence data resulted in a 
matrix with 4676 base pairs (bp), out of which 2676 were 
conserved and 1155 were parsimony-informative. This 
matrix was used to perform all phylogenetic and biogeo- 
graphic analyses and was partitioned by gene and coding 
positions into 19 sections (Table S3). These data were not 
saturated as observed by transitions/transversions rate 
(3.2) compared with genetic distance. 



F. F. Roxo et al. 

A Bayesian analysis was performed with 50,000,000 
generations that resulted in 50,001 trees, of which the first 
5000 were discarded and the remaining 45,001 were used 
to construct the consensus tree (Fig. 2a). The MP analysis 
with the "new technology search" (Goloboff 1996, 1999) 
resulted in a single most parsimonious tree (Fig. 2b) with 
a consistency index = 0.509 and a RI = 0.682. All nodes 
were supported by Bremer indexes with values ranging 
from 1 to 156, and most of the nodes were also sup- 
ported by high bootstrap values. 

Figure 2(a) (ML and BI analyses) and Figure 2(b) 
(MP) illustrate that Neoplecostominae is only monophy- 
letic when including the genus Pseudotocinclus, with 
considerable statistical support (BS = 100 with ML and 
P = 1 with BI, BS = 100 and B = 23 with MP). 

Three lineages are recognized within Neoplecostominae. 
The first lineage (clade A) is composed of Pareiorhina 
rudolphi, P. cf. rudolphi, and Pseudotocinclus (BS = 100 
with ML and P = 1 with BI, BS = 100 and B = 28 with 
MP), and is distributed across both eco-regions, littoral, 
and interior. The second lineage (clade B) is composed of 
Isbrueckerichthys, Pareiorhaphis, Kronichthys, and the 
species Neoplecostomus ribeirensis (BS = 100 with ML and 
P = 1 with BI, BS = 100 and B = 32 with MP) and is 
almost exclusive to littoral. The third lineage (clade C) 
is composed of Pareiorhina carrancas, P. cf. carrancas, 
Pareiorhina sp. 1 a new genus, and all species of the genus 
Neoplecostomus, except Neoplecostomus ribeirensis (BS = 75 
with ML and P = 1 with BI, BS = 51, and B = 5 with MP). 

Relaxed molecular clock and historical 
biogeography 

The mean substitution rate for the Neoplecostominae data- 
set estimated using Beast was 0.808% per MY. The origin 
of Neoplecostominae was 36.9 million years ago (95% 
HPD: 22.3-52.1) with a crown age estimate for the 
Neoplecostomus of 21.3 million years ago (95% HPD: 11.6— 
31.1) (Fig. 3). Within Neoplecostominae, clade A origi- 
nated 29.5 million years ago (95% HPD: 17.4-42.1). The 
ancestral lineage of Pseudotocinclus tietensis reach the inte- 
rior from littoral drainages 8.6 million years ago (95% 
HPD: 2.7-15.4), and Pareiorhina cf. rudolphi 1 9.5 million 
years ago (95% HPD: 3.0-17.0). Results of Lagrange and 
molecular clock analysis suggest that the splitting between 
clade B and C is 26.7 million years ago (95% HPD: 15.4- 
38.1, Fig. 3), and during that time the ancestral of clade C 
reached the interior from littoral drainages (Fig. 3). Within 
clade B more recent lineages also reached the interior from 
littoral drainages: that is the ancestor of Pareiorhaphis ves- 
tigipinnis 3.8 million years ago (95% HPD: 1.0-7.3), the 
ancestor of Isbrueckerichthys sp. 1 4.8 million years ago 
(95% HPD: 2.3-7.8), the ancestor of I cf. calvus 5.9 million 



2442 



© 2012 The Authors. Published by Blackwell Publishing Ltd. 



F. F. Roxo et al 



Biogeography of Catfish Neoplecostominae 



(a) 



re 



N^ifile: y~: 



P b 10:4 3 



Neoplecosiornus sp 7_33421 

Neoplecostomus sp 3_29331 

— Neoplecostomus langeanii_27990 
. Neoplecosiornus Dandeiranle_i36i6 

1 Neopleccstomus sp. 5J3428 
1 — Neoplecosiornus sp. 4_29S32 
L Neopfecoslomus sp. 3_33403 
|0D p Neoplet ostornus paranensisj 7441 
iT-Neoplecostorrus paranensis_17444 
. ^Neoplecostomus sp. 1_32377 
Neoplecostomus so 2 33440 



- Neopleccstomij:»9<vi_9r01 



,-_j4:;jj 



a Neoplecosiornus s; 
1 ^Neoplecosiornus cf. 6otiicatii_6C49 

"i- Neoplecostomus t)0lucatu_34832 
— Neoplecosiornus corurnlja_9*23 

Neo pi et ostomy; iVan: s.coens.i:_3lt 'I I 

Neoptecoslomis mi crops J 593 

Neoplecosiornus espiritosaritensl3_1524: 

Pveiwnnia c arra nca s_3 7560 
Paieioihiriacarr™cas_37S59 
Paieiornina cf. c arra nca s_1 025 B 
Pareiorhlna sp. 1J37565 
Pareiorhina sp. 1_37564 

Newgenus and species 2_37571 

u . — Paieiorhaphis vestigipinriis_3935 
lcpjn — Patetohaplils rystrh<_736"3 

I Pareiorriaptiis stendacrineri_79B9 

I Pareloihapriis spiendens_S257 

- Krorilchthys la 
jnicfitriys subteres_6334 
- Kronii: hlrr/s lieylandij 5096 
nelitlwssp 

lefVeilchlliv" ■ 




(b) 




Neoplecosiornus paianensisj 7441 
Neoplecoslomus paranensis_1 7444 
Neoplecosiornus sr. i_32377 — 

Neoplecoslomus sr. 2_3344Q ™ 
Neoplecosiornus sp. 3.33403 
Neoplecoslomus sp. 4_29B32 _ 
M^ople.: os'oinuE. L^nd'.-ii anteJS61BW 



ll.-up.-.u-.vi.inu;. sp. ij_lu24i; 

Neoplecostomus sp. 7_3M 21 
Neoplecosiornus langeanii_2799G 
I leopleios!omus sp. 3_29991 
Neoplecosiornus sp. 9_10268 
Neoplecosiornus cf 
Neoplecosiornus be 
Neoplecosiornus sel 

N'/uplivos'omuj yapo_970l 9 
Miople.: os'i:muE: coiumb.3_ , 342 3 
Neoplecosiornus rriic rops_7S9? 
Neoplecosiornus espnitosatiterssisj -.Mi 
Neplecostomus franciscoensis_31681 
Parelorriiria carrancas_37559 £ 
Pareiortiiria carrancas__37560 
Pareiortiiria Cf carrancas_10256 




1_37564 
1_37565 

species 2J7570 



P -j : ■■ 'ii ii : 
P -j '« nr. 

Kronichthys lar_erla_83Q4 
Kroriichthy; ?ijl:t*ie:_6^j4 
ivonichthy s lie,- 1 sndMSu ¥1 



3S | — Pareiorhaprns \y itri;_7363 
looPT — Pareioriiapliis wesliBiprinis_S935 

±U I Pareioriiapliis steindacliiieri_7989 

I Pa t off ji*v, soK-jeri fe; - 

BB-Leckp'ictfi\? :' csl.us 'ATI* 
■» 1 ebLDckoicrtivi -,o 1 :53o-: 
" BB-LKke'icwrvs 3j9en 17432 

BB-LKkE'icrryD 0MkTOS_:4E52 
IsbrueckericHhys alipiorils_34853 
Neoplecosiornus ril]eirensis_34S37 



H. 1 -i-.t-i it i- -inrrn-ii 1 =it ■ 1 =: 



Cllojiiiirise 

t_32387 J 

r PRIj ] H',l».(j ■„[.,., -..Ji:;. .;■ 

^HjpffilorriiHe 
— Rineloricarla jaran.uerisis_8369 ]Ln 



pEeu-Mui. r i,lu;-;ielvr.iii_l 
Pseudolocinclu;. iuiiuse_7-:.- 
Pareiorhina cf. mdolphi_241 
P.jrel^iiiiiij ■uij.:l:-.hi_3:775 
go nus Find suei. ti. l_3 
Corumbatala cuestae_12191 
PVpopiopoma lnf*peclatuni_PR12 ~J Hjptpioponerti 




U'liwis 



ik Htigt ge 



_is::e.3 ■ 



['.:!■ r...j:- 



- Hvpoitomui nigrcmacuialus_21 309 

- R rii?li: 1 car s jat.3(iuer'iiis_ : 3::6a ~]b 

- Hernlpsillchlrr/sBoUloJ5363 

- Hemipsilic hires pap ila\iis_1 0241 



]H, p., 



Figure 2. (a) Majority-rule consensus tree obtained in ML and Bl analyses. Numbers above branches are bootstrap values from 1000 bootstrap 
pseudoreplicates obtained from ML analysis. Bootstrap values below 50% {-) are not shown. Numbers below branches are posterior probabilities 
obtained in the Bl analysis, (b) Single most parsimonious tree obtained in the MP analysis (CI = 0.509, Rl = 0.682). Numbers above branches are 
bootstrap values from 1000 bootstrap pseudoreplicates. Numbers below branches are Bremer decay support values. Bootstrap values below 50% 
(-) are not shown. 



years ago (95% HPD: 2.9-9.5), and the ancestor of J. alipio- 
nis 10.6 million years ago (95% HPD: 5.5-16.5). Also, 
lineages from interior reached the littoral drainages 
between 3.5 and 24.0 million years ago: ancestral node of 
New genus and species 2 24.0 million years ago (95% 
HPD: 13.3-34.5), Pareiorhina sp. 1 6.2 million years ago 
(95% HPD: 2.3-11.2), ancestral node of Neoplecostomus sp. 
9 3.5 million years ago (95% HPD: 1.6-5.5), the ancestor of 
Neoplecostomus microps, and N. espiritosantensis 10.2 million 
years ago (95% HPD: 4.7-16.5). Thus, considering our 
Lagrange results (Fig. 3), we observed more lineages of Neo- 
plecostominae reaching the interior from littoral drainages 
(six events) than littoral from interior drainages (four events). 

Discussion 

Headwater captures dispersal hypothesis 

Our analysis suggests that the ancestor of the Neoplecos- 
tominae originated in littoral drainages 36.9 (22.3-52.1) 



million years ago (Fig. 3), based on the results from our 
relaxed uncorrelated molecular clock analysis. The same 
coastal origin was proposed by Chiachio et al. (2008). 
Ribeiro (2006) proposed some interesting biogeographic 
patterns exhibited by the freshwater ichthyofauna on the 
Brazilian crystalline shield and the Atlantic coastal drain- 
ages. The first pattern was defined by sister-group 
relationships between the endemic ichthyofauna of the 
Brazilian coastal drainages and adjacent shield that 
includes intermediate degree of inclusive and both 
sister-clades underwent some radiation. Ribeiro (2006) 
postulated that the members of genera Lignobrycon, 
Rhinelepis, Spintherobolus, and Triportheus, the tribes 
Aspidoradini and Glandulocaudini, and the subfamilies 
Cheirodontinae and Sarcoglanidinae are examples of this 
pattern. This distribution of species is similar to that found 
in our results to subfamily Neoplecostominae, in which 
species of clade B are almost exclusive to littoral drainages and 
those of clade C are almost exclusive to the interior. 
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Figure 3. The red coloration in the figure indicate the interior drainages and green the littoral drainages in southeastern of Brazil. Numbers in 
the map refer collection points listed in the Supplemental Material 1. The tree is the beast chronogram tree from 100 million generations, 
showing diverge age in subfamily Neoplecostominae, calibrated on the origin of Loricariidae {85-90 million years ago) based on Lundberg et al. 
(2007). Scale = millions of years before present. The Maximum-likelihood inference of geographic range evolution in Neoplecostominae, 
performed using Dispersal-Extinction-Cladogenesis implemented in Lagrange v 2.0. L, littoral; I, interior is indicated in the tree. Asterisk indicates 
when the ancestral of clade C reaches the interior between 26.7 (15.4-38.1) million years ago. 



Another pattern of distribution of fish fauna recognized 
by Ribeiro (2006) is related to the families Doradidae and 
Trichomycteridae, in which a taxon of ancient fish 
inhabiting littoral rivers is a sister group to a much more 
inclusive, derived, and widespread taxon. Additionally, all 
species of the subfamily Delturinae (which has been sug- 
gested by Armbruster 2004 as sister group of all other 
loricariids, excluding the genus Lithogenes) are distributed 
in littoral rivers (Rio Paraiba do Sul and Rio Doce). This 
reinforces the hypothesis that the coastal drainages of 
Southeastern Brazil represent an important origin for 
many lineages of Neotropical catfishes. However, this 
cannot be considered a general rule applicable to the ori- 
gin of all Neotropical fishes. Weitzman et al. (1988), 
while working with morphological phylogeny for Glandu- 
locaudini, proposed an opposite evolutionary scenario 
where ancestral species inhabit the Paraguay and Parana 
basins (interior), whereas more derived species are 



distributed along coastal rivers (littoral). However, Weitz- 
man's hypothesis for the Glandulocaudini has yet to be 
tested using additional data and fossil-calibrated analysis. 

Our results also suggested that many species of Neople- 
costominae are distributed throughout elevated regions of 
Southeastern Brazil (Fig. 3). Species of this subfamily are 
known to require waters with high oxygen concentrations 
(Langeani 1990), and it is possible that they have 
difficulty adapting to rivers in lower elevations where the 
oxygen concentration is low. Hoorn (1993) discusses two 
great marine incursions that occurred between 17 and 
20 million years ago and another between 10 and 
12 million years ago. These incursions most likely inun- 
dated lower elevation rivers of southeastern Brazil up to 
the La Plata basin (Hoorn 1993). This suggests that the 
ancestral Neoplecostominae species inhabiting rivers that 
were inundated by those marine incursions are likely to 
have gone extinct; if not entirely, then perhaps at the 
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local population level. Such events may have also pre- 
vented subsequent lineages from successfully recolonizing 
and invading rivers affected by past marine incursions. 

Within the Neoplecostominae, our phylogenetic 
analyses identified three lineages. The first is clade A, 
with species distributed throughout both littoral and 
interior drainages, originating 29.5 (17.4-42.1) million 
years ago. The second is clade B, with species almost 
exclusive to littoral drainages, and the third is clade C, 
which is almost exclusively distributed in interior drain- 
ages. Clades B and C form sister groups and originated 
26.7 (15.4-38.1) million years ago. Albert and Reis 
(2011) reported that several headwater capture events 
occurred between 15 and 28 million years ago among 
the Rio Tiete, Rio Paraiba do Sul, Rio Sao Francisco, 
and Rio Ribeira de Iguape basins. This process is likely 
to have influenced ancestral fish movements throughout 
adjacent drainages, such as the clade C, reached the 
interior from littoral drainages (Fig. 3), and ancestral 
lineages reached the littoral from interior drainages (e.g., 
the ancestor of a new genus and species 2 24.0 
(13.3-34.5) million years ago). Our results also suggest 
that more recent lineages of Neoplecostominae have 
reached the littoral from interior through head water 
captures, such as the ancestor of Pseudotocinclus tietensis 
8.6 (2.7-15.4) million years ago and that of Pareiorhina 
cf. rudolphi 9.5 (3.0-17.0) million years ago. Ribeiro 
(2006) related that the activations of the ancient faults, 
responsible to head waters captures, occurred several 
times in the Serra do Mar formation, an area of south- 
east Brazil where these species are distributed. Riccomini 
and Assumpcao (1999) suggested that the Continental 
Rift of Southeastern Brazil underwent significant restruc- 
turing (upper Miocene, Pliocene, Pleistocene, and 
Holocene), resulting in the reactivation of ancient struc- 
tural lines. This finding reinforces the hypothesis that 
southeast Brazil has undergone intensive geological 
activity and that the activations of ancient faults could 
have resulted in headwater captures between adjacent 
drainages during several periods of its geological history. 

Within clade B, we analyzed five species of the genus 
Isbrueckerichthys, and the results from the Lagrange 
analysis (Fig. 3) suggest that the ancestor of this genus 
inhabited littoral drainages. However, Isbrueckerichthys sp. 
1, I. cf. calvus, and I. alipionis inhabit interior drainages. 
This result suggests that headwater captures could be 
responsible for the current distribution of the species 
inhabiting interior drainages, with ancestral lineages 
reaching the interior from the littoral 4.8 (2.3-7.8), 5.9 
(2.9-9.5), and 10.6 (5.5-16.5) million years ago. Mon- 
toya-Burgos (2003) previously suggested that headwater 
capture events occurring between coastal drainages and 
the Upper Rio Parana (over 4.2 million years ago) have 
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also influenced the current distribution of the genus Hyp- 
ostomus. 

For the genus Pareiorhaphis, we included only four of 
the 19 nominal species (Pereira et al. 2007, 2010) in the 
present study. However, the species not included herein 
occur throughout almost all portions of the Oriental 
Coastal Basin and only a few species are found in interior 
drainages. Specifically, our results suggested that 
P. vestigipinnis could have reached interior drainages 3.8 
(1.0-7.3) million years ago, further supporting the hypoth- 
esis that headwater captures occurred several times in the 
geological evolution of southeastern Brazil. Species of the 
genus Kronichthys are all from littoral drainages, suggest- 
ing that this genus could have originated there. 

Within clade C, our results suggest that Neoplecostomus 
originated in the middle Miocene 21.3 (11.6-31.1) million 
years ago (all species of Neoplecostomus, except 
N. ribeirensis) and is older than the origin of any other 
Neoplecostominae genus (see Fig. 3). Despite the genetic 
diversity found in Neoplecostomus, these species are 
known to be morphologically conserved and difficult to 
distinguish (Langeani 1990; Zawadzki et al. 2008). How- 
ever, it remains unknown whether this morphological 
similarity is due to historical mixture between interior 
and littoral drainages that mixed populations. Our results 
suggest that after the colonization of interior drainages, 
some Neoplecostomus ancestors reached the littoral from 
interior drainages (most likely via headwater captures), 
such as the ancestor of Neoplecostomus sp. 9 around 3.5 
(1.6-5.5) million years ago and the ancestor of Neopleco- 
stomus microps and N. espiritosantensis around 10.2 
(4.7-16.5) million years ago. Furthermore, two lineages 
closely related to the genus Neoplecostomus reached the 
littoral from interior drainages: the ancestral of the new 
genus and species 2 24.0 (13.3-34.5) million years ago 
and the ancestor of the Pareiorhina sp. 1 6.2 (2.3-11.2) 
million years ago. 

Earth climate oscillation dispersal 
hypothesis 

The sea level and tectonic deformations of the continental 
platform have changed drastically in the last 100 million 
years ago, but more intensively in the last 10 million 
years ago as a result of climatic oscillations (Rossetti 
2001; Miller et al. 2005; Miiller et al. 2008; Santos et al. 
2008; Zachos et al. 2008; Albert and Reis 2011). After the 
Middle Miocene Optimum Climate, the climate of the 
Earth changed quickly to a drier and colder resulting in 
the Antarctic Reglaciation (Petit et al. 1999; Zachos et al. 
2008). It was a period of great climate oscillation that was 
intensified throughout the Pliocene, Pleistocene, and 
Holocene (Suguio et al. 1985) resulting in multiple 
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regressions and transgression of marine waters (Hoorn 
1993; Lovejoy et al. 1998; Lundberg 1998). In ages of 
marine regression, animals that lived in the coastal rivers 
might have had contact with adjacent drainages through 
river deltas. These connections would have been impossi- 
ble during periods of high sea level (Montoya-Burgos 
2003; Hoorn and Wesselingh 2010). This hypothesis can 
help to explain the great diversity of species distributed 
throughout the coastal drainages of Southeastern of Brazil 
and specifically the dispersal of Neoplecostominae genera 
through unconnected adjacent coastal drainages. Our 
results suggest that most of the lineages in Neoplecostom- 
inae appeared in the last 10 million years ago (Fig. 3), 
during a period of high variation in the sea level. These 
variations are often related to oscillation of climate in 
Miocene, Pliocene and Pleistocene periods (Ab'Saber 
1979; Colinvaux et al. 2000; Wilkinson et al. 2006; Zachos 
et al. 2008). Such events also affect species living in the 
interior drainages and could have influenced the dispersal 
of Neoplecostominae species within the hydrographical 
basin. Montoya-Burgos (2003) also related this hypothesis 
of delta connections during marine regressions to explain 
dispersal of lineages of genus Hypostomus through 
adjacent coastal drainages in Amazon and northeastern of 
Brazil. 

In conclusion, we hope that the previously discussed 
hypotheses can help in the conservation policies for 
freshwater fishes in southeastern Brazil, an area of vital 
importance in the evolution of Loricariidae and other 
Neotropical fish groups. Phylogenetic and biogeographi- 
cal knowledge continues to be essential for the elucida- 
tion of the magnitude of species richness and species 
evolution and distribution patterns. However, as the 
mechanisms that drive the great diversification of 
Neotropical fish remain understudied and significant 
phylogenetic information is lacking for many groups 
occurring in eastern Brazil, this region still deserves sig- 
nificant research effort. 
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